Skyrmion crystals in the pseudo-spin- 1/2 Bose-Einstein condensates 
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Exact two-dimensional solutions are constructed for the pseudo- spin- 1/2 Bose-Einstein conden- 
sates which are described by the coupled nonlinear Gross-Pitaevskii equations where the intraspecies 
and interspecies coupling constants are assumed to be equal. The equations are decoupled by means 
of re-combinations of the nonlinear terms of the hyperfine states according to the spatial dimensions. 
These stationary solutions form various spin textures which are identified as skyrmion crystals. In 
a special case, the crystal of skyrmion-antiskyrmion pairs is formed in the soliton limit. 
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I. INTRODUCTION 



II. METHOD 



A skyrmion is a particle-like topologically nontrivial 
soliton 1 which is studied in a variety of research fields, 
especially in the condensed matter physics such as the 
quantum Hall effects^O, the liquid crystals^, and the he- 
lical ferromagnetics'. It is noted early that skyrmions 
in three dimensions (3D) have physical properties of 
baryons while skyrmions in two-dimensional (2D) play 
an important role in condensed matter systems. In re- 
cent experiments, it has been demonstrated that the 2D 
skyrmion spontaneously appears as the ground state in 
the helical magnets^ and in the quantum Hall systems at 
fillings slightly away from unity. 

Since the experimental realization of quantized vortices 
in alkali atomic Bose-Einstein condensates (BECs), phe- 
nomena related to the internal degrees of freedom are ar> 
pealing for the investigation of topological objects^^. 
Furthermore, precise spin manipulation techniques have 
been developed to prepare topological spin structures of 
interest and provide unique opportunities to study their 
stability and dynamics**"**!. The skyrmion-skyrmion in- 
teraction may lead to the formation of a Skyrme crys- 
tal. Meanwhile, the periodic skyrmion crystal configura- 
tions could be used to model nuclear matter^M*^. Such 
a state was recently observed in a neutron scattering 
experiment^. 

In this paper we attempt to seek exact stationary solu- 
tions to the 2D coupled nonlinear Gross-Pitaevskii equa- 
tions (GPEs) which describe the pseudo-spin- 1/2 BECs. 
The intra- and inter-species interactions are assumed to 
be equal and the system is exposed to a uniform exter- 
nal field. We propose a method to decouple the GPEs 
by re-combinations of the nonlinearity and the spatial 
dimensions. Periodic solutions are obtained which are 
shown to form the skyrmion crystals in the pseudo-spin 
representation. In Sec. II we describe the method. The 
main results are displayed in Sec. III. Section IV contains 
a brief summary. 



We consider the 2D pseudo-spin- 1/2 BEC in a uniform 
external potential (V(r) = 0) by assuming the intra- and 
inter-species interaction strength 711 = 712 = 7 and the 
two-component atoms have the same mass mi = 1122 = 
m. The Hamiltonian then has the "pseudospin" SU(2) 
Symmetrel By adopting the units of h = m = 1, the 
stationary GPEs for the mean-field order parameter \I/ = 
(^i,^2) T are 

/^;(r) = [~V 2 + 7(1^1 (r)| 2 + |^ 2 (r)| 2 )]^(r), (1) 

where jij (j = 1,2) denote the chemical potentials. In 
order to solve the coupled Eqs.(TJ), we decompose the 
wave functions according to the dimensions as 



ip 1 {x,y)=X 1 {x)Y{y) 
Mx,y) = X 2 (x)Y(y) 



■iX(x)Y 1 (y) 
■iX(x)Y 2 (y) 



(2) 



where Xi(x),X2(x),X(x) and Y\ (?/), ¥2(2/)? Y(y) are real 
functions of a single variable. By imposing the restric- 
tions of 



\Yi(y)\ 2 - 



\X 2 (x)\ 2 
\Y2(y)\ 2 



(3) 



where B\ and B2 are real constants, the total density can 
be seperated as 



p(r) = B{X\x) + BiY\y). 



(4) 
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Substituting Eq.Q into Eq.Q and equaling the real and 
the imaginary parts, respectively, we obtain a set of ordi- 
nary differential equations for X(x) and Y(x) which can 
be solved self-consistently by making use of the unique 
properties of the Jacobian elliptical functions*^. The de- 
tails are described in the next section. 

According to the relations of three non-singular Ja- 
cobian elliptical functions, sn, cn, and dn, we construct 
three forms of solutions to the Eqs.Q. They exhibit peri- 
odic crystals of the skyrmions or skyrmion-anti-skyrmion 
pairs, respectively. The periods of the wavefunctions are 
k = 4 • jK(m), where K(m) is the complete elliptic inte- 
gral of the first kind with modulus < m < 1 and j is a 
integer which denotes the number of periods. 
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III. STATIONARY SOLUTIONS 
A. Type A 

We first consider the following form of solution, 

J ip\(x, y) = Asn(fcx, m)cn(%, m) + iBcn(kx, m)sn(%, m) 
1 ^2(^7 2/) = Dcn(kx, m)cn(ky, m). 

(5) 

This wavefunction automatically satisfies the periodic 
boundary conditions in the x- and ^-direction due to the 
properties of the Jacobian elliptical functions. By sub- 
stituting Eq.([5| into the stationary Eq.Q and using the 
identities between the Jacobian elliptic functions, we ob- 
tain the relations: 



A 2 = B 2 



1 



D 2 



and 



Ml = I/c 2 (2-3m 2 ) 
H2 = k 2 (l-2m 2 ), 



7 



km 
"^42" 



(6) 



(7) 



(8) 



From 7 < 0, we note that this solution applies to the BEC 
with attractive interactions. The difference between the 
chemical potentials \±i (i — 1,2) of the two components 
can be induced by a Zeeman energy. The total number 
of atoms is related to the amplitude of the wavefunction 

by 



2A 2 



1 



E(m) 
K(m)'' 



(9) 



where E(m) is the complete elliptic integral of the second 
kind. 

Figure 1(a) and (b) display the checkerboard distri- 
butions of density of each hyperfine state for A = 1, 
7 = -11.3670, /ii = 28.4175, /i 2 = 22.7340 and m 0.5. 
The phase of ipi shown in Fig. 1(c) reveals a periodic 
vortex- anti- vortex pair structure. The phase diffe rence 
implies a particle current density which is defined b} * 18 * 29 * 



1 

2% 



^[^(v^)-(v^*)^ 



(10) 



j=l,2 



Fig. 1(d) displays the vortex- anti- vortex structure of the 
total particle current. 

In order to explore the topological structure of the 
state, we examine the normalized spinor x( r ) which is 
defined by \I/(r) = y p(r)x(r). The spin vector S for the 
condensates are S(r) = x^X? with a the Pauli matrices. 
Fig. 2 (a) displays the spin texture which forms a skyrmion 
crystal. The topological charge or the Pontryagian index 
of the skyrmion is an invariant!^, 



Q = q(r)dxdy, 

J cell 



(ii) 
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FIG. 1: (Color online) The density profiles of the ipi (a) and 
ip2 (b) hyperfine states in the state |5]). The parameters are 
A = 1, 7 = -11.3670, /ii = 28.4175, fi 2 = 22.7340 and 
m = 0.5. The two-species have different particle number and 
different chemical potential, (c) The phase distribution of 
the complex ipi component, (d) The total particle current 
density. The number of period of the wavefunction j = 1. 



where the topological charge density 



q(r) 



47T 



S- (d X S X dyS). 



(12) 



Fig. 2(b) shows distribution of q(r) which exhibits a pe- 
riodic structure. We note that in one period of the state 
([5]), there are four skyrmions. Each unit cell has a total 
topological charge Q = 1. By the way, the anti-skyrmion 
crystal can be obtained by simply take the complex con- 
jugacy of the state 
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FIG. 2: (Color online) (a) The spin texture S(r) of the state 
forms a skyrmion crystal. The number of period j = 2. 
(b) The distribution of the topological charge density q(r). 
The total topological charge of each cell is unity. 
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FIG . 3: (Color online) The same as in Fig.l for the state 
(13). The parameters are A — 1, 7 = 11.3670, /ii = 56.8351, 
/i2 — 51.1516 and m — 0.5. 



and 



A 2 = B 2 = -D 2 , 
2 ' 



/ii = ±/c 2 (2 + ra 2 ) 



/i 2 = /c 2 (l + ra 2 ) 



2^ j 



(14) 



(15) 
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B. Type B 

Next we consider the solution for 7 > which is ful- and 
filled by the following configuration, 

J ip\(x, y) = Asn(fcx, m)cn(%, m) + iBcn(kx, ra)sn(%, m) 
1 ^(z^y) = Dsn(kx,m)siti(ky,rn). 

(13) 

We obtain the relations: 



FIG. 4: (Color online) The same as in Fig. 2 for the repulsive 
case of the state (13). The total topological charge in each 
cell is Q = — 1. 



C. Type C 

We finally show an example of solution that is generally 
not a skyrmion as m < 1 while tends to be a skyrmion 
in the soliton limit (m — >• 1). The solution is of the form 

J i/)i(x, y) = Asn(kx, ra)dn(%, ra) + iBdn(kx, ra)sn(%, m) 
1 ^2(^7 2/) = Ddn(kx, m)dn(%, ra). 

(18) 

We obtain the relations: 

1 



A A =B Z 



-m z D z 



Mi = 



1 1,2 



k 2 (2m 2 - 3) 
~fc 2 (m 2 -2), 



7 = To 



*1 

A 2 ' 



(19) 



(20) 



(21) 



Obviously, the solution applies to the attractive BECs. 
The total number of atoms is related to the amplitude of 
the wavefunction by 

>E(m) 



N = 2A 2 



K(m)' 



(22) 



7 



k 2 m 2 



Figure 5 display the results as those in Fig.l for the 
state (18). The parameters are chosen as A = 1, 7 = 
-45.4681, 111 = -56.8351, jn 2 = -7U5691 and m = 0.5. 



(16) Despite of the similarity, the state (18) is essentially dis- 



The total number of atoms is related to the amplitude of 
the wavefunction by 



N= 2A^ EM 
ra 2 K(rn) 



(17) 



Figure 3 display the results as those as in Fig.l for 
j = 1. The parameters are A = 1, 7 = 11.3670, \i\ = 
56.8351, H2 = 51.1516 and m = 0.5. The spin texture 
forms a periodic crystal of anti-skyrmions as shown in 
Fig. 4(a). Our calculation shows that each anti-skyrmion 
has a total topological charge Q = — 1 (Fig. 4(b)). 



tinct to former two states. To clarify this point, we show 
in Fig. 6 the spin textures for various modulus ra. We 
find that the total topological charge defined by ( pT| ) for 
each domain is not an integer. In a period, the spin tex- 
ture does not form a skyrmion since the spin vector S(r) 
does not rotate a complete circle around the z-axis, as 
shown in Fig. 6 (a) and (b). Nevertheless, in the solitonic 
limit we find the spin texture still form a skyrmion-anti- 
skyrmion crystal in comparison to the former solutions. 
Fig. 6(d) reveals that the total topological charge in each 
period is alternatively Q = ±1 as m 1. This configu- 
ration is shown in Fig. 7, where the charge density forms 
a dipole crystal. 
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IV. SUMMARY 
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FIG. 5: (Color online) The same as in Fig.l for the state ( 18 ) 



The parameters are A = 1, 7 = —45.4681, jii — —56.8351, 
/i2 — —79.5691 and m — 0.5. 



(a) m=0.1 




(c) m=0.999, 



In summary, we have presented a class of exact solu- 
tions to the coupled GPEs which describe the pseudo- 
spin-1/2 BECs. The spin textures exhibit periodic crys- 
tals consisting of skyrmions (Q = 1 or Q = —1) or 
skyrmion-anti-skyrmion pairs. 
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FIG. 7: (Color online) The same as in Fig. 2 for the state ( 18 ) 
in the limit of m — > 1. The total topological charge in each 
cell Q = ±1 which forms a dipole crystal. 



FIG. 6: (Color online) (a)-(c) The spin texture of the state 
(T8j for m = 0.1,0.8, and the solitonic limit m = 0.999. (d) 
the total charge defined by (11) versus the modulus m. As 



m — > 1 the spin texture becomes a skyrmion-anti-skyrmion 
crystal. 
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